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Abstract 

Charge is transported through superconducting SSS single-electron transis- 
tors at finite bias voltages by a combination of coherent Cooper-pair tunneling 
and quasiparticle tunneling. At low transport voltages the effect of an "odd" 
quasiparticle in the island leads to a 2e-periodic dependence of the current 
on the gate charge. We evaluate the I — V characteristic in the framework 
of a model which accounts for these effects as well as for the influence of the 
electromagnetic environment. The good agreement between our model cal- 
culation and experimental results demonstrates the importance of coherent 

Cooper-pair tunneling and parity effects. 
PACS numbers: 73.40. Gk,74.50.+r,73.40.Rw 
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The single-electron tunneling (SET) transistor has proven an ideal system to display the 
concepts of "single electronics". In this device an island is coupled via tunnel junctions to 
the leads. The island potential can be modulated by a capacitively coupled gate voltage Vg. 
In transistors with normal- conducting islands and leads the current depends e-periodically 
on the gate charge Qo = CgVg. Recently much attention has been devoted to single-electron 
transistors with superconducting islands and normal leads (NSN) and entirely superconduct- 
ing systems (SSS). In a superconducting island, where Cooper pairs form the condensate, 
the addition of one extra electron - the "odd" one - costs the gap energy A [|l|. Hence the 
physical properties of the system depend on the parity of the charge number in the island, 
and the I — V characteristics is expected to be 2e-periodic in the gate charge. The first clear 
signature of 2e-periodicity was observed by Tuominen et al. 0|. However, no satisfactory 
explanation of these experimental results has been provided until now. In the meantime 
a number of experiments both in NSN transistors and in SSS transistors f^-^ have 
demonstrated a rich variety of phenomena and show good agreement with theoretical re- 
sults. It is now also well understood ^ that the difficulties in observing the 2e-periodicity 
arise from the extreme sensitivity of the even-odd difference to effects of the electromagnetic 
environment which create non-equilibrium quasiparticles. 

Parallel to the experiments the theoretical description of systems exhibiting parity effects 
made rapid progress. In Ref. the authors present an equilibrium model which accounts 
for the temperature dependence of the even-odd asymmetry. A kinetic model was developed 
in Ref. [ffOl to describe transport. In this framework the I — V characteristics of NSN 



transistors can be derived For SSS transistors, there exists a well-developed theory 



for the current-biased system [12,5,13|. The voltage-biased system has been considered in 



the absence of parity effects while in Ref. ||^ only resonant Cooper-pair tunneling has 
been studied. In this paper we investigate a model which includes Cooper-pair and single- 
particle tunneling as well as parity effects. First we study the I — V characteristic of an SSS 
transistor in the absence of external impedances. In the transport voltage range cV rijE(j 
even-odd effects are observed {Ec denotes the scale of the charging energy, see below). We 



then discuss examples for the relevant transport processes. In order to compare with the 
experiment, we account for the influence of the electromagnetic environment. This rather 
complex model explains the experimental results of Ref . . 

We consider a SET transistor (see Fig. 1) with superconducting electrodes and island 
(with energy gaps A) below the crossover temperature T^r = A/kBlnNeff where parity 
effects can be observed. Here Nejf = 2Nj{0)\/2nAkBT is the effective number of states 
available for the odd quasiparticle and A''/(0) is the density of states (per spin) in the 
island. For the moment we ignore the effect of the external impedance. 

The system can be described by a sum of Hamiltonians for the left and right electrode and 
the island, the charging energy and the tunneling Hamiltonian H = Hi+Hji+Hj+Hch+HT- 
We will treat the Josephson tunneling non-perturbatively. Hence we start from the model 
Hamiltonian 
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and we will account for the quasiparticle tunneling in perturbation theory. Here Q = 
{ni — nr)e denotes the island charge, and we defined a total charge which has passed through 
the system Q = (rii + nr)e (with respect to a reference state). Further C = Ci + Cr + Cg is 
the total island capacity and Qq = CgVg the offset charge (here we assume Ci = Cr). The 
first two terms describe the charging energy and the energy gain in tunneling due to the 
transport voltage. They are diagonal in the basis of charge states \Q,Q)- The typical scale 
of the charging energy is Ec = /2C. The last term describing Cooper-pair tunneling with 
coupling energy Ej is off-diagonal and, thus, mixes the charge states \Q,Q)- Therefore, the 
eigenstates of Hq are linear combinations of these states 

l^fe) = E<ml'^e,me) . (2) 

n,m 

A dc-current requires a dissipation mechanism. In the case of zero external impedance 
the quasiparticle tunneling can cause transitions between different eigenstates \^k)- It is 



accounted for by 

K = E T^^2|Q + e,Q + e)(Q,g|cIc, + h.c. 
+ E tI:;}\Q - e,Q + e){Q,Q\cl,Ck + h.c. 

k£l,p'GR 

The tunneling matrix elements in the left/right junction T^j!^^ are related to the conductance 
of the junctions by l/Ri/r = 47re^/^ Nl/r{0)Ni{0)\T'^''^''^\'^, where we approximated T^,!''^ ^ 
jf junction resistances are large compared to the quantum resistance Ri/r > Rr = 
h/e"^ the transition rates can be calculated by the golden rule 

S-^oToVi-expl-^i/ABT) ; 

■ E E |(M/;|Q±e,g±e)(Q,Q|vI/,)|2 . 
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Here l'^^ is the well-known / — ^ characteristic for quasiparticle tunneling between 
superconductors and Sij is the energy difference between initial and final state. In order to 
allow for the parity effect, we have to include the escape rate -^esc of an odd quasiparticle in 
the island ||T0| 

lesc — 
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In order to determine the dc-current, we follow the procedure described in Ref. Il^ 



First, we determine the eigenstates of Hq either in perturbation theory (as we shall dis- 
cuss below) or numerically taking into account a sufficient number of charge states. This 
procedure converges for not too large Josephson coupling energies Ej < Eq- Given the 
eigenstates \^k) we calculate the rates in Eq. (^, which then enter a master equation 
dtPk = J2n=ik{Pn^n^k — Pk^k-^n) foT the probabilities Pk to find the system in the k-th eigen- 
state. The stationary solution dtPk = is sufficient to evaluate the dc-current 

^ = ^ E PkTk^nii^nlQl'^n) - ■ (6) 

k,n^k 



The results are shown in Fig. 2(a). We used the parameters A = l.SEc, Ej = 0.17Ec, 
Ri/r = -R ~ Rk, lesc = 2.5 ■ 10^^(i?C)~\ which correspond to those in Ref. 0. Here 
E^j denotes the Ambegaokar-Baratoff expression for the Josephson energy. In Eq. (|I|) the 



generahzed Josephson couphng energy in the presence of charging effects enters, which 
in the present case is larger than E^j by roughly 20 %. We note that the I — V characteristic 
is 2e-periodic and observe a rich structure deep in the subgap region. For transport voltages 
eV^2.bEc the 2e-periodic features disappear and the current becomes e-periodic in Qq 
again. This is not surprising since on a current scale / ^ e7esc the unpaired quasiparticle 
in the island looses its importance. 

The basis of eigenstates of Hq is "appropriate" for the problem. On inspecting our nu- 
merical procedure we find that for low transport voltages only a few (two or three) states 
Y^k) are noticeably populated. Similar behavior is found in systems without coherent tun- 
neling like NNN or NSN transistors. Therefore, we can calculate transition rates and the 
current analytically if we know the eigenvalues of Hq and the corresponding eigenstates. To 
this end, we determine the coefficients in Eq. (^ by using perturbation theory in Ej. 
Away from certain resonant situations, the fc-th eigenstate has only one coefficient a^^ of or- 
der unity, whereas all other coefficients are considerably smaller. To fix ideas, let us consider 
the state |\l/o) in the range of gate charges Qq G [0,e/2]. In this eigenstate the most likely 
charge state is \Q = Oe, Q = Oe), i.e. ~ 1. Due to coherent tunneling of one Cooper 
pair, there is a small amplitude 0^.3 ^2 °^ ^j/ to find the system in the charge states 
\Q = ±2e, Q = ±2e). Since also several Cooper pairs can tunnel coherently, the system can 
be, e.g., in the charge state |2e, 6e) with a very small amplitude a^g {Ej/EcY. 

At resonance lines, however, it is possible that this amplitude is much larger. Let us 
consider the solid straight line in Fig. 2 (b), which is given by 

3eV = 4Ec{l-Qo) • (7) 

Along this line the charge states |0e, Oe) and |2e, 6e) have the same energy, i.e., three Cooper 
pairs tunnel resonantly there. The resonance results in a drastically increased amplitude 



cx {Ej/Ec). 

A transition from to another eigenstate can occur if it is energetically favorable and 
the matrix element of the final state with |\I^o) according to Eq. (|^) is nonzero. On analyzing 
which transitions due to quasiparticle tunneling are energetically favorable, we find that a 
process 

l^o) ~ |0e, Oe) — > l^i) ^ |le, 7e) (process a) 

is possible. Out of resonance the rate of process a) is of the order (Ej/Ec)^. In a narrow 
strip (whose width is characterized by Ej) around the resonance line Eq. (0), however, we 
find 

' Ej/2 Y (Ej-^ 



The line in Fig. 2(b) corresponding to Eq. (^ marks the most significant resonance in the 
I ~V characteristic. We are, thus, led to the conclusion that the dominant transport process 
in the subgap region is tunneling of quasiparticles accompanied by simultaneous tunneling 
of several Cooper pairs. Due to this combination enough energy is gained to overcome the 
quasiparticle tunneling gap 2A. The importance of this type of transport mechanism was 



first noted by Fulton et al. [jT^- Although the rates for these processes in general are small, 
they are considerably enhanced in situations, where resonant transfer of Cooper pairs is 
possible. 

So far we have studied the conditions for the system to leave the initial state. However, 
a dc charge transport through the system requires cycles, after which the island returns to 
a state equivalent to the initial one. The simplest version is a two-step cycle of subsequent 
transitions of the same type in the left and right junction. Such cycles dominate in NNN 
or NSN transistors at low bias voltages. The cycles which lead to the pronounced features 
in Fig. 2 arise due to two-step cycles as well, but the second step is different from the first 
one. The transition completing the cycle which starts with process a) is 

l^i) ^ |le, 7e) — > 1^2) ~ |0e, 12e) (process b) , 



i.e., quasiparticle transfer accompanied by only two Cooper pairs tunneling. The latter 
process is not in resonance and, therefore, the rate is Ff, oc (Ej/Ec)^. Whereas off-resonance 
the process a) is the bottleneck for the current, at resonance the process b) has the smaller 
rate. This explains that at the resonance the current increases by roughly two orders of 
magnitude. 

Another interesting feature in the I — V characteristic is the shoulder-like structure 
between the high resonances for gate charges Qo E [e/2,e]. It is directly related to the 
escape rate ■jesc of the odd particle. The first step in the relevant cycle is a process similar 
to process b) 

l^o) ~ |0e,0e) — > 1^3) ^ |le,5e) 

with a rate oc (Ej/Eq)'^, which is relatively large (as discussed before). The current, however, 
is limited by the second step 

1^3) ^ |le,5e) — > 1^4) ^ |0e,6e) . 

This is a pure quasiparticle transition without Cooper pairs, which can occur because the 
escape rate jesc has no gap. In Fig. 2(b) it is seen that the cycle sets in for transport voltages 
eV > AEciQo — 1/2) (the dashed straight line). This is exactly the condition for the odd 
quasiparticle to gain energy on leaving the island. 

So far we have considered the ideal case of a vanishing external impedance Z. In order 
to compare our results with experimental data it is necessary to account for the effect of the 
electromagnetic environment. An external impedance gives rise to incoherent Cooper-pair 
transitions. For not too low transport voltages the rate of these transitions is given by (see, 



e.g., Refs. 0,113) 

rg^4l(.,|Q/2|..)r ^_^^);^^^ . (9) 

For appropriate parameters they lead to a pronounced resonance structure in the I — V 
characteristic pO| . 



Furthermore, in an experiment the temperature of the environment is not necessarily the 
same as the electron temperature [pT|,^. High temperature fluctuations in the environment 



induce quasiparticle tunneling and, thus, can cause qualitative changes in the I — V char- 
acteristic. In order to take into account this effect, we have to add to the single-electron 
tunneling rate Tj^j- the term 

for Eif < 2 A. The Ohmic resistance Re describes the environment, which fluctuates at the 
high temperature T^^Ec- Since there is no detailed information about the electromagnetic 
environment in Ref. 0, we flx the parameters by comparing with similar experiments. 
In Fig. 3 the I{Qo) dependence for several transport voltages is plotted. The calculation 
reproduces remarkably well both the general shape of the experimental curves and the order 
of magnitude of the current. We have cut the current resonances for the bias voltages 
eV = 167 fiV, 200fiV, since the maximum differential conductance, which can be observed 
in an experiment, is limited p]. In a model calculation like ours we cannot expect perfect 
correspondence between theory and experiment. The reason is that the system is very 
sensitive to even small amounts of quasiparticle tunneling which can overcome a threshold 
for cycles in lower order of Ej/Ec- (We have discussed a similar effect in connection with 
the shoulder-like structure). Different relevant temperatures in the environment (which is 
very likely) can cause subtle changes in the I — V characteristics. In this case Eq. (|10]) is 
only a rough approximation. Therefore, we have to adjust the parameter for different 
bias voltages. Finally we mention that one can speculate about processes of higher order in 
the quasiparticle tunneling, e.g., co-tunneling or coherent two-electron tunneling across one 
junction, accompanied by Cooper pairs. However, we come to the conclusion that in the 
present case processes of higher order in the quasiparticle tunneling can be neglected. There 
are two possibilities: i) The process creates excitations. Then the rate has a gap > 2A and 
becomes important only for higher transport voltages, where also flrst order quasiparticle 
tunneling is possible, ii) The process does not create excitations. But then the phase space 
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is reduced by at least one small factor 1/(A'"/(0)A) [|l|. 

Summarizing, we can say that our non-perturbative model adequately describes charge 
transport in voltage biased SET transistors under various conditions in a large range of 
transport voltages. Due to the parity effect the I — V characteristics are 2e-periodic in 
the gate charge for not too large transport voltages. The dominant transport mechanism 
is quasiparticle tunneling, accompanied by coherent tunneling processes of several Cooper 
pairs. An external impedance at low temperatures causes incoherent Cooper-pair transitions, 
whereas high temperature fluctuations in the environment "poison" the ideal structure, thus 
restoring e-periodicity of the current. 

We acknowledge stimulating discussions with R. Bauernschmitt and A. Rosch and thank 
P. Joyez for a copy of his thesis. This work is supported by "Sonderforschungsbereich 195" 
of the Deutsche Forschungsgemeinschaft. 
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FIGURES 

FIG. 1. The SSS transistor, n^/^ is the number of electrons which have tunneled through the 
left and the right junction. 

FIG. 2. a) I— y characteristic of a SSS transistor for vanishing external impedance (parameters 
see text), b) Contour plot of the same data. There are three dashed contour lines in the current 
range 7 = 0... 10~^ e/{RC) and 20 lines for I < 10~^ e/{RC). A pronounced resonance is found 
along the straight solid Hne 3eV = 4:Ec{l — Qo). The straight dashed line eV = AEc{Qo - 1/2) 
marks the edge of the shoulder-like structure. 

FIG. 3. The current I{Qo) through an SSS SET transistor, including both resonant Cooper-pair 
transitions due to an external impedance Z and high temperature fluctuations in the electromag- 
netic environment. Parameters are: Ec = ISO/ieV, A = 240/ieV, Ej = 30/icV, Ri + Rr = bOkfl, 
jesc - 2 • 10S-\ ReZ(0)=80J), Re = O.Sn (for V < 123/xV), Re = Sn (for V > 150/xV), re=3.6K. 
We have assumed a slight asymmetry of the junction resistances Ri : R^ = 4 : 3. 
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